NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



FINAL SCIENTIFIC REPORT 


Organization: 

The Regents of the University of California 
University of California, San Diego 
La Jolla, CA 92093 

Department: 

Department of Mathematics 

Principal 

Investigator: 

M. Rosenblatt 

Title: 

Three-dimensional bl spectra and the mean 
rate of transfer of energy 

Agency and 
Grant No. : 

National Aeronautics and Space Agency 
NSG-2376 

Period: 

3/1/79 - 10/30/80 


(NASA-CR-164217) ES1IMA1ING Nal-2*a22 

THREE- DEHLNSICNAL ENERGY 1RAKSEER IN 
1SOT BOPIC TURBULENCE Fina* Scientific 

Report, 1 Mar. 1979 - JO Cct. I960 Uncias 

(California Uuiv., Sax. Die-jo, La JOi.ia.) G3/o7 4212 7 


Estimating three-dimensional energy transfer in isotropic turbulence 


K. S. Lii 

Department of Statistics 
University of California, Riverside, CA 92502 

K. N. Helland 

Institute for Pure and Applied Physical Sciences and 
Department of Applied Mechanics and Engineering Sciences 
University of California, San Diego, La Jolla, CA 92093 


M. Rosenblatt 
Department of Mathematics 

University of California, San Diego, La Jolla, CA 92093 


1. Introduction 


The research reported on in this paper is a continuation of that consid- 
ered in previous work (Lii et al. (197 6), Helland et al. (1979)). Here the 
sf. .mpt is to get a three-dimensional spectral view of grid-generated turbulence. 
This is done by setting an x-wire probe at a fixed position, and setting two 
single wire probes at a number of locations in the same plane perpendicular to 
the mean flow in the wind tunnel. The locations of the single wire probes are 
determined by pseudo-random numbers (Monte Carlo). Second-order 
spectra and cross -spectra are estimated. The assumption of isotropy 
relative to second-order spectra is examined. Third-order spectra are 
also estimated corresponding to the positions specified. A Monte Carlo 
Fourier transformation of the downstream bispectra corresponding to 
integration across the plane perpendicular to the flow is carried out 
assuming isotropy. Further integration is carried out over spherical 
energy shells. The object is to get an estimate of the spectral transfer 
function that indicates the rate of decay of energy. 
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2. Basic theoretical background 
2. 1 Energy transfer equations 

Let v(x, t) = (Vj.V 2 .Vj) be the velocity field as a function of loca- 
tion x » (Xj.x^. Xj) and time t. In the model of homogeneous turbulence. 
v(x, t) is regarded as a random solution of the Navier-Stokes equation 


gxlEJl + v . 

at 


7 v = 


1 . 2 - 
— 7p + V 7 V 
P 


and the continuity equation 


7 • v = 0 


that is stationary in x. Here 7 is the gradient operator, p is the density, 
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p is the thermodynamic pressure, v is the kinematic viscosity, and 7 
is the Laplacian. Given finiteness of second order moments, it follows that 
v(x, t) has a Fourier representation 

v(x, t) = /e ik x dz(k,t) , 

Here k = (kj.k 2 .kj) is the wavenumber vector. Existence of appropriate 
moments is assumed whenever required. The spectral matrix f(k, t) is 
given by 

f(k,t) 6^- p dk = E fdz(k,t) dz*(k',t)} 

= (E {dz (k, t) dz Q (k', t)} ; a, P = 1, 2, 3 ) . 
a P 

where the asterick is used to denote the transpose conjugate. From now 
on the time dependence will be implicit and the explicit appearance of t 
will ho suppressed. The third-order momento 
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Q a g(k,k') dk dk' = Im E f k a dz a (k-k') dZg(-k) dyk')} ; 


oi 1 3 = 1 , 2 t 3, 


The energy associated with wave number vector k is f..(k) (repeated 
indices are used to imply summation over i = 1, 2, 3). 

The equation 

Jr(z f ii (i '>) = /Q(k.k'>dk'-v|k| 2 f..(k) 

relates second- and third-order spectra and describes the rate of change of 
energy with time. 

Q(k.k') = y q Q <k,k') 
a.p a ’ 8 

is considered to be the net mean rate of energy transfer from dk' to dk. 

The continuity equation implies that 

Q(k.k') = -Q(k',k). 

2. 2 Symmetries 

At times we shall refer to the three spatial locations as x , x x " . 

At other times, using the homogeneity, we shall refer to the first point as 
(0, 0, 0), the second point as x and the third point as x / . There are a 
number of symmetries which hold under the assumption of isotropy 
for third -order moments or third-order spectra. We shall first 
derive these for third-order moments. In Fig. 1 the coordinate axes are 
indicated as well as velocities at three spatial locations. At location x the 
velocity is given in the x^ and x^ directions. At locations x' and x" 
only velocity in the x^ direction is given. The velocities indicated in Fig. 1 
correspond to the experimentally measured velocity components. First 
notice that the three-point correlation: 
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C ul <x.x\ x ') * Efv^x) v^x') v^x*)} 


= C {«;.* , x \' x \ » x 2* x 2' x 2' x 2 1 x 3‘ x 3* x 3' x 3 ^ 


by stationarity in the spatial variables (homogeneity). Isotropy implies 
that we can change the coordinate system (x^.x^. x^) by permuting any of the 
coordinate axes. Interchanging the x^ and x^ axes we obtain 

E £ v 2 (x 2 ’ X r X 3 ) v 2 (x', x ;.x') v 2 (xj. xj.x')] 

3 C 222^ x 2' x 2‘ x 2’ x 2 ; x l’ X r X 1* X 1 ; X 3* X 3’ X 3 _X 3^ * 

Thus 


C lll <X l' X l* X l'V X 2' X 2‘ x 2* x 2 ’ x 3* x 3' x 3’ x 3 ) 


C 222 (x 2‘ x 2' X 2" X 2 ; x I“ x i' VV* x 3 _x 3’ x 3' x 3 ) 


C 333 (X 3' X 3* VV x 2‘ x 2' x 2 _x 2 : VV 


by interchanging x^ and x^ axes. Also 
C 2ll^ x ’ X> ' * ^ 

= E {v 2 (x l( x 2 , x^) v (Xj, x 2 ', x j) v ^ (x* » x", x") } 

= C 2ll (x !' X r X r X i ; X 2' X 2* X 2" X 2’ X 3‘ X 3’ VV 
= C 122 (X 2 _X 2' x 2* x 2 ’ x r x r x r x l ; x 3' x 3 ’ x 3‘ X 3 ) 


(by interchanging x^ and x 2 axes) 
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* C 3, *;-V x 3 -x 3 ’ x 2' x 2’ x 2- x 2' 

(by interchanging and axes) 

8 ^233^ X 3’ X 3' X 3*" £ 3’ X 2~ X 2* x 2 _X 2’ X l" X l’ X l’ X l^ 
(by interchanging x^ and x^ axes) 

= C 133^ x 2' x 2’ y ?" x 2 ; X 3 -X 3’ X 3 _X 3 ' X l‘ X l’ X 1* X 1^ 
(by the permutation (x^.x^.x^) to (Xj,Xj,X2)) 


= C 322 (x 3' x 3' VV X [ _X 1* x l‘ x l : x 2* x 2‘ x 2 _X 2 ) 


(by the permutation (x^.x^.x^) to (x^, x^.Xj)) 


Now, let 


B agp (k,k / ) dkdk' = E f dz^k-k') dzp(-k) dz g (k') 3 


A simple heuristic argument implies 


dz (k) ^ 

a 


1 ( - ik • 

(2tt ) 3 i 


v (x) dx 

a 


so that 


B a8P ,,t - C '> 


f -i(k-k') • x f ik* x' , 

/e v a (x) dx • / e Vp(x')dx 

/ -ik'*x ff 

e v fl (x")dx" 

1 f ik* r - ik'* r ' 

ur 




C aB8 (r ’ r ' )d?d? ' 


where r = x x , r ' = x" - x. 


The symmetry relations among imply 


B 111 (k,k')*B 111 (k 1 ,k'; k 2 * k 2 #; k 3 ' k 3 '> 


B 222^ k 2* k 2 5 k j* k j # * k 3* k 3 # ) 


= B 3 33 ( k 3» k 3 i k 2 * k 2 5 *V k l # ) 


and 


B 2n°V k i' 1 


k 2’ k 2 


k 3 ', = 

B l22 (k 2' 

k *• 

V 

k *• 

V 

k 3‘ 

k 3» 

s 

B 311<*V 

k '• 

k 3 . 

V 

k 2 . 

k 2> 

= 

B 133 (k 2* 

k '• 

V 

k 3 . 

V 


k l» 

= 

B 233 (k 3* 

V 

k 2' 

k' • 
k 2 ’ 

k l* 

V 

= 

B 322 (k 3 

. k^ ; 

k l 

k '• 
’ k l ’ 

k 2 . 

k 2> 
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3. Monte Carlo and Computational Methods 

Velocities ^(0,0,0; t), v 2 (0, 0, 0; t), v^O.x^Xj*, t) and v^O.x^.x'; t) 
were measured as functions of time t at the fixed locations (0,0,0), (0, x 2 , x^h 
and (O.x'.x'h The x^ and x^ axes represent downstream and perpendicular 
directions, respectively. The velocities as functions of time t can be 
interpreted aB functions of spatial coordinate x^ by using Taylor's 
hypothesis x^ = Ut, where U is the mean downstream velocity. This is 
essentially a "frozen-flow" approximation where the velocity is measured 
at a fixed location, and the turbulent velocity field is swept past at a 
uniform speed U. The frequency spectrum of the fluctuating velocity can 
be interpreted as the wave -number spectrum of the fluctuating velocity 
component in the downstream direction (Batchelor, 1953). The 
choice of (x 2 » x^) and (x 2 , x') is determined by a Monte Carlo method 
such that each coordinate x^. x^, x^ , x^ is uniformly distributed on 
±1,0 inches. There are 100 such pairs ((x^. x^), (x^.x^)). Later, this 
was augmented with seven more pairs with each coordinate uniformly 
distributed within ± 0. 333 inches. The seven pairs were added to enhance 
the resolution at small spatial separations. 

Let 


b age (k r k i : x 2’ x 3 : x 2 * x 3 ) 


C f i(k x - k* x / ) 

— ~"*2 ll e ^ 11 E ( v a (0. 0.0) VgUj.x^Xj) Vg(xj,x 2 . Xj)} dXj dxj 

(2rr) JJ 
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and 


N- 1 


F (k,x) 
a 


x) * A v (x , x_, x t) exp { -itk 2 tt/N) 
Jri*. a i <• * 


t = 0 


F a (k l 5 x 2* x 3 ) 8 F a (Uk l ; X 2’ X 3 ) 


where U is the mean downstream velocity. Then 


b aPf3 (k r k ' ; x 2‘ x 3' X 2* X 3 ) 


F^ikj-k^'; 0,0) F^-kjj x^, Xj) Fp(kj # ; x^.xp 


Suppose we have N pairs (x^\x^), j = 1,..., N obtained by 

Monte Carlo. Then 


B app (k,k,) = E f dz a (k ' - k,) dz g<- k > dz 0 (k')} 


“T / e 

(2tt) 


ik - r - ik'r' 




~b f 1(k 2 r 2 + k 3 r 3 ' k 2 r 2 - k s' r 3 n 

(2rr) 

• [jfjTexp {«!,,», • kj'r,')) C aB9 <V r ,'s ; r y rpd.j dr,'] 


- dr 2 dr 3 dr 2' dr 3 


*77 /j® rexpfl(k 2 r 2 * k 3 r 3 • k 2 r 2 " k 3 r 3 ,) 
(2tt) 


■ b aPe (k l’ k l' ; r 2' r 3 : X V r 3 ,dr 2 dr 3 dr 2 dr 3 


9 


We estimate B oa (k, k 7 ) by 

app 7 


W ^'* 95 


/*\ jl. y b (k k /. x u) x (j). x o) x a), 

V 3 / N(I) ^ aSn’ 1 * 2 ,x 3 ’ x 2 ,x 3 ; 

• exp{i(k 2 x^ + k 3 x^ -k'x'^L k^ x 7 *^) } 


ife-(2/3i 4 y b . k x <j) x o). r ,o) ,o). 

N(E) aPP'V 1 ’ 2 ,X 3 * X 2 ’ x 3 ’ 


j € E 


• exp (i ^ 2*2 ^ + k 3 x 3 ^ * k 7 x 7 ^ ‘ k'x 7 ^^)} 


3 3 


for a = 1, 2 and 8 = 1 where N(I) = 11 = number of pairs 
<x«\x«>s x 2 ' ( ^,xf) which lie in 1 



E = l-l, 1 ] 4 - 1 , N(E) = 96 . 


Because additional points were taken in 1, different weights are applied here. 
The details of the Monte Carlo locations are given in Sec. 5. 

If k =(kj,k 2 . k^), k 7 = (k| , k^ , k ^ then by the methods described 
above, estimates are obtained of 

B m (k,k') = B ul (k l ,kJ; k 2 ,k 7 ; k 3 , k 7 ) 

and 

B 2n (k,k 7 ) = B 2u (k l ,k l 7 ; k 2 .k 2 7 ; k y k 7 ) 
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By the discussion in See. 2 it follows that 


E 22l (k,k') * 1 1^2‘ ^2 * k j* k /5 k 3 * k 3> 

B 333 (k. k') = B nl ( k 3 * k 3 k j» k jj '• kj»kj) 

Bj 22 (k,k') * B 2ll^ k 2' k 2* k l ,k {* k 3» k 3^ 
B 3n (k, k') = B 2n (k l .kJ ; k 3 ,k' ; k^kp 

Bj 33 (k, k') * B 2u (k 2 , k' ; kj, k' ; kj.kp 

B 233^ k * k ^ * D 21 l^ k 3‘ k 3 # ' k 2’ k 2 ’ *V k i ^ 
B 322 (k, k') * B 2u (k 3’ k 3 ; k i ,k i #: k 2 ,k 2 ) 


Now 


Q(k 


k.k') = X Q a 8 (k ' k,) 


a , 8 


can be determined from 


Q (k. k') = k Imli (k. k') 

a r ci a p B 


Estimates of bias and variance of B ou (k, k') are given in the appendix. 

a p p 

To simplify the final presentation and further reduce the variance 


a Monte Carlo shell integration was carried out. We compute 


with k^ , k'^ independent and uniformly distributed on spheres with 
radii r and r' ( respectively. 

The net rate of spectral energy transfer at | k | = r is 

P( | k | ) = P(r) = / P(r, r ') dr' 
m 

* £ P(r, r/) . 

i - 1 

where the r ' are given by an auxilliary Monte Carlo procedure. This is 
to simulate the uniform integration on the sphere. 


4. Numerical Implementation 

The single wire sensor locations (0, x^i x^) and (0, x'. x' ) were 
determined by a Monte Carlo technique. The IMSL subroutine GGUW for 
generating uniformly distributed random numbers was used to generate 
4 x 100 points in the square ± 1 in. by ± l in. and 4x7 points in a square 
±1/3 in. by ±1/3 in. Velocities were measured at locations ((0,0,0), (0, x^, x^)» 
(O.x'^.x'^)) for j = 1, . . . , 107 and are shown in Fig. 2. Eleven of them have 
all coordinates of length less than 1/3 inch. An x-wire located at (0, 0, 0) 

* t 

was used to measure both downstream and perpendicular velocities , v^ 
and single wires located (0, x^, x^) and (O.x^.x^) only downstream velocities 
vj.vj 3 were measured. Compute 


b au"v k r 4 x r x 2 J - x 3 J > * b in<V k i'i 


= WV 1 °’< 0) F 1 (-k J ; x^.xj) Fj(kj # ; x.,' j ,x 3 ' j ) , a = 1.2 


where 


N- 1 


F ^(k;x,, x ) = y v (0, x,, x ; t) exp { -itk 2rr/N } 

cx ^ ^ 


with N = 4096 for 


k = 26,76 976 

k j = - 97 6, . . . , - 26, 26, 7 6, . . . , 97 6 


1 3 




T 


We obtained 


~ b a,,<W * 


kj+25 k|f 25 

1 V V . j 

50x50 ^ b all (m,n) 

m = k^ - 25 n = kj - 25 


for tx * 1 p 2. 9 j = 1 $ • . • » 1 07 as our basic computational data base, 


Im D m (k.k ') 
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^ u)j Im bj^ikj.kp • cos {k^x^ + k^xj - k^' - k^ x^ } 


is computed where 


N«* 


if j C I 


f N(K) lf •> € E * 

In this computation only the imaginary part of ^^(k *k*) is considered. 


The real part of bj is zero by the isotropic assumption and this is 
confirmed by the statistical resolution of the data. Similarly, 


Im iv,<k. k') 

- I I 

107 


^ tu . Re b\ (k . k ') • sin [k,x^ t k x j - k * x '■* - k / k'/* ] 
•‘T'j J 21111 2 2 33 22 33 


since lm h^ = 0 under the isotropic assumption. All other im R (k, k') 

1 1 ex H 

computations are carried out in a similar fashion using the symmetry condi- 
tions indicated in Sec. 3 giving 
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To do t tic spherical shell integration, we use the GGSP11 subroutine in the 
International Mathematical and Statistical Library to generate uniform points 

k^, k , j = 1 2000 on spheres with radii r and r'. 

Finally, 


Q(r.r') = Q( | k | , | k'|) 


10 




10 " m 

m 1 


where 


Q (r, r ') = -r 
m 200 


1 00 m 
\' 

j = lO0Tni-l) + 1 


Q(k J , k /j ) 


represents estimated net energy transfer from the spherical shell with 
radius r ; to spherical shell r . This format was chosen so as to permit 

m 

estimating the variance of Q(r,r ). 


5. Experimental Methoda 

The experimental conditions were chosen to match the conditions used 
previously at the University of California at San Diego by Helland ct al. (1979). 
These conditions are nearly identical to a number of studies carried out by the 
authors and also by a number of other investigators. The experiment was 
performed in the 7 6 cm by 7 6 cm by 9 m test section of the low-turbulence 
wind tunnel in the Department of Applied Mechanics and Engineering Sciences. 

A biplane grid of round, polished dural rods was located 2. 4 m from the end 
of the contraction section. The grid had a mesh size M of 5. 08 cm with rods 
0.953 cm diameter. A DISA model 55P1 x-wire sensor measuring two com- 
ponents of velocity (v and v ? ) was oriented with the plane of the x parallel to 
the wind tunnel walls and was located at a downstream distance X/M = 48, 
where X is the longitudinal distance from the plane of the grid. Two single 
hot-wire sensors DISA model 55P01 were mounted so that each sensor could 
be moved with two degrees of freedom to any point in the transverse plane 
intersecting the x-wire sensor. 

The x-wire sensor was held fixed on the tunnel centerline by a sym- 
metrical airfoil extending horizontally through a slot in the side wail of the 
tunnel. The x-wire sensor was attached to a shaft passing through the center 
of the airfoil. The shaft permitted rotating the x-wire in the vertical plane 
thereby providing a method for calibrating the x-wire sensor over a range of 
yaw angles. A large protractor mounted outside the tunnel was used to read 
the yaw angle to about ± 0. 1 degree of arc. One of the single wire sensors was 
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mounted on a second airfoil support mounted outside the tunnel. A vertical 
slot in the wind tunnel wall permitted vertical motion of the airfoil, while the 
horizontal positioning was obtained by moving the airfoil back and forth 
through the tunnel wall on a set of bearings. The vertical slot was sealed 
by a sliding plate system which prevented air flow out of the wind tunnel. 

The other single wire sensor was mounted on a standard machinist's height 
gauge which was modified to hold the hot-wire support and moved in the 
vertical transverse direction. The height gauge was mounted on a Unislide 
base which moved in the horizontal transverse direction. 

The position of each single wire sensor was measured relative to the 
location of the x-wire sensor location x = (0,0,0) by four dial indicators 
which could be read to ±0.0001 inches. The dial indicators were adjusted 
to be direct reading in the following manner: First the two single wires 
were moved so they were in the same longitudinal plane perpendicular to the 
wind tunnel floor passing through the plane of the x-wire sensor. This was 
accomplished by observing the four wires through a Gaertnor traveling 
microscope. The two horizontal dial indicators were then set to read zero. 
Next the vertical positions of the two single wire sensors were measured 
with the traveling microscope, and the two vertical dial indicators were set 
to read the measured values. The single wire locations were measvired 
before and after each series of experiments. The repeatability was good to 
about ±0.001 inch. Tighter tolerances in positioning the hot-wires would 
require a more rigid support system. 
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A pitot tube pressure sensor operated by a MKS Baratron electronic 
manometer was used to read the mean speed U. The 5 (im diameter by 
1. 25 mm sensing length hot-wires were operated by four DISA 55M01 
constant-temperature anemometers, and the anemometer outputs were 
linearized by four DISA 55D10 linearizers. A Honeywell 7 600 FM analog 
tape recorder was operated at 3-3/4 in. /sec for the hot-wire calibrations 
and at 15 in. / sec for the turbulence data giving an effective bandwith of 5 kHz. 

Calibrations were performed on the hot wire sensors before and 
after each experiment with the grid removed from the wind tunnel. During 
preliminary testing the x-wire probe was set to give symmetrical yaw 
response about the wind tunnel centerline. Each hot-wire was linearized and 
the sensitivities were adjusted to be nearly identical. This enhanced the 
usual sum and difference decomposition approximation for x-wire calibration, 
but this was not an essential step because the analog sum and difference 
signals were not used to separate the anemometer voltage signals into the 
longitudinal (v^) and vertical (v^) velocity components. Instead the x-wire 
sensor was calibrated by setting the mean wind speed using the pitot- Baratron 
instrument to one of [6. 0,6. 5,7. 0,7. 5, 8. 0,8. 5,9.0} m/sec and then 
rotating the x-wire sensor through angles [ 0, ± 1 , ± 3, ± 5, ±7 } degrees. 

The wind tunnel speed was then changed, and the x-wire sensor repeatedly 
yawed through the set of angles. Approximately 63 pairs af mean speeds 
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and yaw angles were used to define the full two-dimensional calibration of 
the x-wire. This calibration technique minimized the assumptions needed to 
interpret the output of x-wire sensors and we were not required to assume 
that the sensors were identical or mounted exactly at 90 degrees to each 
other. The calibration data were analyzed by standard linear regression 
techniques to obtain the constants (a, b, c's) for the calibration equations to 
first-order in anemometer voltage 

Vj (0, 0, 0; t) = a + b e^t) + c e 2 (t) 

v 2 (0,0,0; t) = a' + b'e^t) - c'e^t) 

v^O, x 2> x 3 ; t) = a" + b" e^t) 

V i (O’ x 2 * Xjl t) = a " + b"e 4 (t) 

where e. is the linearized hot-wire voltage of hot-wire channel i. 

Channels 1 and 2 correspond to the two x-wire voltages, and channels 3 and 
4 corresponded to the two single wire voltages. The assumption of first- 
order dependence in the calibration equations was tested by performing a 
regression which included the second-order terms involving e. and compar- 
ing the results to the first-order analysis. An additional check was made by 
comparing some of the bispectra calculated from an identical portion of the 
data from a first-order calibration to one done using the second -order 
calibrations. The differences were found to be negligible in both cases. 

The sets of uniformly distributed points shown in Fig. 2 were grouped 
in such a manner that the experiments could be performed efficiently. For 
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example, the coordinates were grouped so that the sensor locations falling 
in the first quadrant were referred to as Case 1, and sensor locations fall* 
ing in the first and second quadrants were referred to as Case 2, and so on. 
Each pair of coordinates was assigned a number within each case such as 
Case 1. 1, 1.2, etc. A check was performed on the computer generated 
coordinates which assured us that the distribution of points was uniform 
within the standard deviation expected for the sample size. 

The turbulence data for the 107 pairs of single sensor locations 
were measured in a series of seven experiments. The data were taken at 
night between 6 PM and 7 AM to ensure that the wind tunnel temperature 
variation did not exceed ±0.1 °C. Data were recorded on analog tape for 
a ten minute interval at each pair of sensor locations. The entire series of 
experiments was set up and performed during a 6 week period during the 
Fall of 1978. 

The hot-wire calibrations for all seven experiments were digitized 
onto a single digital tape using an Interdata 7/1 6 minicomputer based data 
acquisition system having 1 2-bit resolution. The turbulence data were 
digitized at a sample rate of 4098 samples/sec, 4096 samples /channel/ 
record and 250 records per sensor location. Analog low-pass filtering at 
a cut-off frequency of 2 kHz was used during the digital recording process 
to minimize aliasing. A standard format was adopted putting two sensor 
locations on each digital tape giving a total of 54 digital tapes required to 
contain the turbulence data. The digital tapes were shipped to the 
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Lawrence Berkeley Laboratory CDC 7 600 computer located in Berkeley, 
California, and processed remotely from San Diego using both remote 
dial-up access and a remote batch terminal located in the NORPAX building 
at the University of California, San Diego. Most of the computer output, 
including graphical computer output, was written to standard 48X power 
microfiche. The microfiche made it possible to process and store the 
several thousand pages of printed and graphical computer output. 
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6. Discussion 


One -dimensional Spectra 

Power spectra, probability densities, and third- and fourth-order 
moments were computed for the 107 pairs of hot-wire positions in order 
to check the quality of the data. Several positions have been selected 
and sample results will be shown only for those locations. Table 1 
specifies the selected locations, and Table 2 shows some typical 
moments and characteristic length scales for the flow condition? 

The power spectra are shown in Fig. 3 for the three locations. The 
power spectra for channels 1, 3, 4, measuring the velocity in the 
x direction, are almost indistinguishable from each other as ex- 
pected. This is a result of the homogeneity of the flow field in the trans- 
verse direction. There is no variation in the spectral shapes among the 
three locations which represent very close, intermediate, and distant 
probe locations. This result suggests that no probe interference is 
apparent in the second-order spectral data. Channel 2 is the spectrum of 
the velocity fluctuation in the x^ direction and should not correspond to 
the spectra in the x^ direction even under the assumption of full isotropy. 
Again there is no evidence for probe interference problems in the x^- 
component spectra. Note that graphs and tables are expressed in frequency 
while the initial derivations of the paper are in terms wavenumber. 

The probability density functions are shown for the selected three 
locations in Fig. 4. The density functions have each been normalized by the 
standard deviation for each channel. All density functions plotted in this 
manner for each location are coincident and nearly Gaussian in shape. 
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Note that again there is no apparent variation in the density curves across 
the three probe locations. While the deviation from Gaussianity appears 
slight, it is important that there be some deviation so that energy transfer 
among fluctuations of various scales may exist consistent with a three- 
dimensional turbulent field. 

An elaborate check on second-order isotropy was made using about 
20% of the hot-wire locations covering the range from |r| small to 
| r| large. From Batchelor (1953) it is easily shown that the cross- 
spectra between the x-wire sensor and a sensor placed at r = (0, r^, r^) 
can be computed from an estimate of the downstream and transverse correla- 
tions f and g 


u f (r ) = E Vj(x) Vj(x + r) 


u 2 g(r) = E v 2 (x) v 2 (x + r) 

where f and g are functions of r = |r|, the separation between the tw 
sensors and u 2 = var(v). Then by isotropy we can write 


and 


R..(r) = 
ij 




2/ f - g 


' r i r j 


+ g 6. 


«) 



rf' 


f ' 


af 

Or • 


Thus the cross spectrum between downstream velocity components 
v, (0,0,0) and v,(0,x-,,x,) can be expressed in terms of f and g 

1 1 w J 
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2 f , , ik,r, 2 r 

. .. . u / If 2 11, u I 

h , (k, , r,, r,) * -y- I - ■= — r, e dr + -r- I 

1, 1 1 2 3 2 tt J 2 r 1 2 tt I 


ik r 

g e 1 1 dr, 


Similarly the cross spectrum between velocity components v^tO, 0,0) and 
v^O.x^x^) is given by 


»2,l (k l’ r 2' r 3> 


_ 00 

3 — [ — 

2nJ r . 


ik l r l 

r l r 2 e dr l 


The correlations f and g are shown in Fig. 5 for the Cases 4. 10 
in both linear and st tnilogarithmic coordinates. The agreement among the 
correlations is excellent with Case 4. 10, and 4. 10 was used to perform 
the isotropy checks which follow. 

The directly estimated cross spectra are compared with the cross 
spectra computed using the assumption of isotropy are shown in r. series 
of figures for the three cases of small, intermediate, and wide separation. 
The results are shown on linear scales in Figs. 6 to 8 to show the com- 
parison at low frequency. The results are repeated on semi-logarithmic 
scales in Figs. 9 to 1 1 where the cross spectra have been multiplied by the 
frequency to show the comparison at high frequency. The agreement is 
generally good, with some significant deviations occurring at the higher 
frequencies. While the assumption of isotropy (and Taylor's hypothesis) 
for these computations appears reasonably well satisfied, it is impossible 
to know the significance of deviations on the integrity of the third-order 
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energy transfer computations where the assumption of full Isotropy is 
required. 

The bispectra b and b. for the three locations are shown 

ill “I 1 

in Figs. 12 to 14. The bispectra do vary with probe locations now since 
they involve three velocities, one from each point (0,0, 0), (O.x^.x^), and 
(0, x^.Xj). The magnitude of each bispectral type does not change much 
with position (all bispectra figures have the same scales), but the "width" 
of the peak does change rapidly with increasing probe separation. The 
larger the separation the smaller the extent of the main lobe of the bi- 
spectral peaks. The variation with probe separation is more readily 
observed in diagrams of the statistical resolution of the bispectra b^^ 

and b_ shown in Tables 3 to 5. In the resolution diagrams, the numbers 
£ 1 1 

represent the rounded ratio of the bispectral estimate to the standard devla 
tion of the estimate. The sign is the sign of the bispectrum. 

The results of the computations for Q(k, k*) are shown in 
Tables 6 and 7 and in Figs. IS- 19. Table 6 shows the magnitude Q and 
the standard deviation of Q i6 shown in parentheses. The truncated ratio 
of the magnitude of Q with its standard deviation is shown in Table 7. 
Figures 15-19 show representative cross-section plots of the data in 
Table 6. The curve s represent energy transfer from a frequency magni- 
tude, f' , to a selected frequency f. For (' < f we expect most of the 
energy transfer to be positive, that is a loss of energy from the low 
frequencies to the higher frequencies; for f # > f we expect most of the 
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energy transfer to be negative, that is, a gain of energy in the higher fre- 
quencies from the lower frequencies. The continuity equation requires 
that Q>0 when f * f ' . This requirement is approximately satisfied in 
Figs. 17 and 18. In Figs. 15, 16, and 19 the Q differs significantly from 
zero when f = f' . This may be due to too broad a smoothing in the spectral 
estimation. In Figs. 17 and 18, the two curves are consistent with the 
expected direction for the energy transfer. These data are in qualitative 
agreement with some calculations of Deissler (1977) who used an expansion 
technique to estimate the shape of the net energy transfer between wave- 
number shells. These computations are in the central spectral range where 
one would expect to have greater reliability because isotropy would probably 
be satisfied in this range (as contrasted with low wavenumber) and we have 
rcsolvability (as contrasted with high wavenumber). 

The physical significance of this investigation is associated with the 
notion of energy transfer among the eddies or scales in the Fourier tran*- 
form domain. Such a transfer of energy in homogeneous turbulence is an 
aspect of the nonlinearity in turbulence. There has been a concept of local 
transfer of energy employed in a number of discussions of turbulence. 
Though the notion has usually not been formalized, it appears to be a claim 
that appreciable transfer takes place only to neighboring wave numbers. 
Since there are three wave number vectors (with sum equal to zero) involved 
in the equation for dissipation and transfer of energy, this means that two 
of the wave number vectors are close to each other and the third is close 
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to zero. The usual assumption has been that the transfer of energy is one- 
sided, from lower wave numbers to higher wave numbers (in terms of 




wave number vector magnitude). This does not seem to be literally con- 
firmed by the results one has that appear to have resolution and that seem 
to be consistent with the usual symmetry Q(| k| , | k'| ) = -Q(| k' | , | k| ), 
results obtained in an intermediate wave number range. It is true that 
one appears to have low variance in the low wave number range but there 
are fears of nontrivial bias here. In the high wave number range one has 
to oe resigned to high relative variance. From the results it seems as if 
there is nontrivial interaction between wave number magnitudes differing 
by a factor of three. 

There are a variety of questions that arise if one were to consider 
carrying out an experiment like this again. It would seem to be more 
appropriate to use two x-wires and one single wire rather than one x-wire 
and two single wires in the experiment. This would make the experiment 
a bit more difficult. However, this would allow one to estimate a 2, 2, 1 
bispectrum as well as the 1,1,2 bispectrum. These two bispectra are 
related theoretically under the assumption of isotropy. However, our 
estimates of the 1,1,2 bispectrum generally had little resolution while it 
is clear in terms of previous experimental results that we would have good 
resolution for the 2, 2, 1 bispectral estimate (Helland et al. (1979)), at 
least for separations close to zero. We could have one fixed x-wire and 
let the x-wire and single wire probe float in the fixed perpendicular plane. 
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If one could get adequate analytic representations of the 1,1,1 and 2, 2, 1 
bispectra as they depend on spatial separation, suggested by the estimates, 
the Q function would be obtained by a direct Fourier transform. 

The Taylor hypothesis has been used to translate time separation 
into the spatial x^ separation. Tests of the Taylor hypothesis with 
separation in the direction have been carried out (Favre et al, (1955), 
Stegen and Van Atta (1970)). In our application we do have separation in 
the perpendicular direction. 

Should the Monte Carlo locations be laid out uniformly over a fixed 
range or set out in a stratified manner? If we have some knowledge about 
the general shape of the function to be integrated, we can use the importance 
sampling scheme to choose a convenient probability distribution which 
approximates a normalized positive integrand. If we lay out the Monte 
Carlo locations according to such a distribution we can get a substantial 
variance reduction (Klcijnen (197 4), Niederreiter (1978)). Also we could 
take advantage of the assumed isotropy and concentrate all locations in a 
subregion using antithetic variates to reduce variance (Kleijncn (1974)). 
Quasi-random points could be used instead of the pseudo- random numbers 
to improve the overall convergent rate in the Fourier transform if the total 
variation of the integrand is not too large (Niederreiter (1978)). 

One of the basic questions is that of better resolution. There is 
no difficulty in appreciably larger number of records (a factor of 10). A 
real difficulty is that of the increased computational cost. With the use of 
modern ultrahigh speed array processes and signal display devices 


2S 


operated by modest laboratory computers, v'he computing cost can be fixed 
and the computing speed is reasonable. 
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Appendix*. Bias and Variance 


In this section the bias and the variance using a Monte Carlo layout 
of the experiment is compared with a deterministic layout. Let 


A a8p (K 


__ f i ( k . x ,+ k . 

<,K') = /e 


X 3 ‘ k 2 X 2 * k 3 X V 


b a?8 (k r k I : X 2* X 3 : X 2’ X 3 ) dx 2 dx 3 dx 2 dx 3 


where K = (k^»k,,k^) , K ' = (kj.k'.k') 


b a8p (k r k V X 2’ X 3 : x 2 ,x 3* 


(2rt) 


i ".' 1 


E {v (0, 0, 0) v (u , x,, xj v Q (uf, x', x' )} du.du' 
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which is a one dimensional bispectrum. An estimate of b „ is given bv 

a ep * 1 
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where the s.'s are such that 
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3 

I , . 0 

i 1 1 
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with 


d (N »,u.) . V X («,/ iai * 

a j J tM> “j 


y. ... 


such that 2* ci . = 0 (mod 2rr). Here 
j = l J 


W N (U 1’ U 2' U 3 ) = “T W 
h 


h H\ 

i h * h * h / 


w is a symmetric weight function such that 


/ 


w(a l* a 2* a 3 ) 6 


x„U, 

. . 1 / 1 

v i = 1 / 


du 2 duj = 1 


It is known that (Brillinger and Rosenblatt (1967)) the bias 


b N = ^ b aBfl^ k l' k P ’ b a6P^ k l’ k 0 


_0 

Ok 


1 


‘WV’'!’' /of U 2 I d »! d 'V 


1 
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0 2 2 

+ — * — j 6 u (k , kj ) f u w(u,u)dudu. 
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+ o(h 2 ) 


if b^^ is twice continuously differ ertiable and appropriate conditions are 
imposed on the weight function w. Under suitable conditions the variance 
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V * r<£ ase (k r l '! >) '77- ‘a 

h N 


(kj > ^( k j) ^( k j + k') / w 2 (uj,u 2 )dujdu. 


°(hW) 


where f is the power spectrum of component a. 

To approximate A ^ • suppose a deterministic uniform grid on 

[-a,a]^, say m^ points, x. . = (x.,x.; x',x'), are given. Then a 

1* Ji K, I i j K x 

natural approximation of A ^ is given by 


m 
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\ m I 


if one can assume twice continuously differentiable. Now 

we assume the independence of measurements at different locations. 
This is reasonable since the measurements at different locations are taken 
at different times. Then 
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Now suppose a uniform random sample of M points from [-a, a] are 

taken, say, x. = (x~., x ; x', x ' ), j = Then a Monte Carlo 

J 2 J 3 J 2j 3j 

approximation of A (K, K') is given by 
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Therefore it appears that when 1 /m is substantially less than h then the bias of 

deterministic layout is better than the bias of Monte Carlo layout, whereas the 

variances are comparable if the number of sample points are the same. When h 

is substantially less than 1 /m then the bias of a Monte Carlo is better than that of 

a deterministic layout and the variances might be comparable. This depends on the 

choice of other parameters. Now the details of the variance will be computed to the 

first order. However, it must be kept in mind that each choice of a set of locations 

corresponds to a complete set of experimental measurements at these locations. 

Further the greater the number of locations, the greater the length and cost of com- 
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putaiions. So even though it is easy to speak of 1 /m small, this amounts to m 
locations ami the experimental and computational cost become enormous. This is 
the great advantage of the Monte Carlo procedure so as to keep experimental and 
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2 

Here O (k ,k') is the variance of Im f ...(k ,kp for all j's. This vari- 

1 1 l t L ” 1 1 

x 

ance is independent of the location because the variance is proportional to the 
product of the power spectra. The previous expression equals 


A f D - C 

where A, B, C correspond to the 3 terms in the preceding equation. 
To evaluate A, drop the subscript j. Then 


a ^ 

J cos (^*2 + ^3 X 3 ‘ ^2 X 2* dx 2 dx 3 dx'dx' 
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Now j cos( a + p - y - 6 ) da d 8 dy d5 
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= j tcos(a+B)cos(Y+6)+sin(a+P) sin( y + fi)] da d B dy 
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J [(cos a cos B - sin a sin B )(cos y cos 5 - sin y sin 5 ) 
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+ (sin a cos 8 + cos a sin B)(sin y cos 6 + cos y sin 6 )] da d B dy d5 
a 

j" (cos a cos B cos y cos 5) da dB dy dfi 


. a 

lince J sin wx dx = 0 for J- 1 * 0. We obtain 
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The values of D and C cannot be determined theoretically since we don't 


know theoretical values of f ...(k.k,) though we have estimated values. 
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Now B and C may be estimated by using the estimates of f 
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Note that B and C are independent of (theoretically). 


In Summary: 
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where this time the imaginary part of b is zero acain by the 
assumption of isotropy. 
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This equals 
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where we use the same notation as before and ignore possible correlations 
among (k-.k^ ; L>^) for different j's. Carrying out the same kind of 
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argument we have; 
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Again A' can be computed from the estimated power spectra which are 
fairly reliable. B* , C* can only be estimated empirically. 
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Table 2. Physical parameters for grid-generated turbulent flow. 


U (m/sec) 

7.77 

V {m/sec) 

0. 002 

M (cm) 

5.08 

d (cm) 

0.953 

r m 

25300 

x/M 

48 

eW 2 1 1/2 /u 

0.0183 

Etv 2 ] I/2 /U 

0. 0142 

Etu 2 ] 1/2 /Elv 2 ] 1/2 

1. 15 


E[u 3 ]/Elu 2 ] 3/2 

0. 0087 

E[v 3 ]/E[v 2 1 3/2 

-0.0186 

Elu^/Elu 2 ! 2 

2.912 

E[v 4 ]/E[v 2 ] 2 

2. 905 

R x 

34.6 

e (cm7sec^) 

787 

ti (cm) 

0.047 

(cm/sec) 

3. 32 

* k = u/2 "\ 

2630 


4 1 


Table 3. Re solvability of one -dimensional bispectra 








































































































Table 5. Resolvability of one -dimensional bispectra 



























































Table 6. Three-dimensional energy transfer spectrum Q(f, f f ) representing the net energy transfer 

from spherical shell of radius f ' to the spherical shell of radius f. The standard deviations 
of the estimate Q is shown in parentheses. 
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Normalized probability density at selected locations g 










Comparison of experimental and derived cross spectra for Case 6. 10. □ , Experimental reel part; 

A , experimental imaginary part; , theoretical spectrum derived from autocorrelations of 

Fig. 5 using isotropy. 


















(,X*J) T Z 0 





Comparison of experimental and theoretical cross spectra for Case 8. 4. □ , Experimental real part 

A , experimental imaginary part; theoretical spectrum derived from autocorrelations of 

Fig. 5 using isotropy. 
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Fig. 15. Net three dimensional energy transfer from 






Fig. 17. Net three dimensional energy transfer from f' to f * 300 Hx 




dimensional energy transfer from f 1 to f * 350 Hz 








